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Abstract: As is known, the so-called Dirac A'-operator commutes with the Dirac Hamiltonian for 
arbitrary central potential V{r). Therefore the spectrum is degenerate with respect to two signs 
of its eigenvalues. This degeneracy may be described by some operator, which anticommutes with 
K. If this operator commutes with the Dirac Hamiltonian at the same time, then it establishes 
new symmetry, which is Witten's supersymmetry. We construct the general anticommuting with 
K operator, which under the requirement of this symmetry unambiguously select the Coulomb 
potential. In this particular case our operator coincides with that, introduced by Johnson and 
Lippmann many years ago. 
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Supersymmetry (SUSY) of the Hydrogen atom is 
rather old and weU-studied problem. We have in mind 
the usual cases: SUSY in non-relativistic quantum me- 
chanics and inclusion of spin degrees of freedom by Pauli 
method as well. In this last case the projection of well- 
known Laplace-Runge-Lenz vector onto the electron spin 
direction plays the role of superchargei. 

Relatively less is known about the Dirac electron, 
although the so-called radial SUSY was demonstrated 
a long-time ago^. As for 3-dimensional case, it was 
shown^i^i^, that the supercharge operator is the one, in- 
troduced by Johnson and Lippmann in 1950 in the form 
of a brief abstract^. As regards to the more detailed 
derivation, to the best our knowledge, is not published in 
scientific literature. Moreover as far as commutativity of 
the Johnson-Lippmann operator with the Dirac Hamil- 
tonian is concerned, it is usually mentioned that it can 
be proved by "rather tedious calculations."— 

The main aim of our paper is a derivation of the 
Johnson-Lippmann operator in a simple and transparent 
manner and simultaneous demonstration of its commu- 
tativity with the Dirac Hamiltonian. 

Below we show that among all central potentials the 
Coulomb potential is a distinguished one. The additional 
symmetry takes place only for this potential. Then we 
show that the operator responsible for that symmetry 
reduces to the Johnson-Lippmann one in case of Coulomb 
potential. 

So, let us consider the Dirac Hamiltonian for arbitrary 
central potential, V^(r): 



H ^ d ■ p 13m -(- V{r) 



(1) 



In this form V{r) is a fourth component of a Lorentz 
4-vector. We mention this fact here because the pure 
Lorentz -scalar potential is also often considered^. 

Let us introduce the so-called Dirac operator—, 



K 



S • ? + 1 



(2) 



where I is the angular momentum vector, a and /? are the 
usual Dirac matrices and S is the electron spin matrix 



E = 



a 
a 



It is easy to show that 

[K,H] 







(3) 



(4) 



for arbitrary central potential, V{r). 

Therefore the spectrum of Dirac equation is degenerate 
with respect to eigenvalues k of the K -operator. As a 
rule, this is a degeneracy with regard to the signs of n, 

We can find an operator A, which could connect these 
two signs. Naturally, such an operator should anticom- 
mute with K, 



{A,K} ^AK + KA^O 



(5) 



If at the same time this operator commutes with 
Hamiltonian H, it'll generate the symmetry of the Dirac 
equation. 

Therefore, we are looking for an operator A with the 
following properties 



[A,H] = 0, {K,A} = 



(6) 



After that we will be able to construct supercharges as 
foUows^i^ 



Qi^A, 



Q2 



AK 



Then it is obvious that 



{Qi,Q2} = 



(7) 



(8) 



and we can construct Witten's superalgebra, where Q\ = 
= h IS a. Witten's Hamiltonian. Now our goal is 
a construction of the A operator. For this purpose at 
first we generalize one theoremi*ii, known from Pauli 
equation to the case of Dirac equation. For the Dirac 
case this theorem may be formulated as follows: 
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Theorem: Suppose V be a vector with respect to the 
angular momentum I, i.e. 

[h,Vj] = iSijhVk 

or, equivalently, in the vector product form one has 

fxV + Vx 1 = 2iV 
Suppose also that this vector is perpendicular to I 

(t- v'^ = ^V-lj = 0. 

Then K anticommutes with operator{Yi ■ V) , which 
is scalar with respect to the total J momentum, i.e. it 
commutes with J = / + 1/2S. 

The proof of this theorem is almost trivial - it is suf- 
ficient to consider the product (S • I) in this and in 
reversed orders and make use of definition of K. Then it 
follows that 

{iC, S • y } = (e • + (e • if = (9) 

It is evident that the class of operators anticommuting 
with K (so-called, K-odd operators) is not restricted by 
these operators only. Any operator of the form 0(S • I) , 
where O commutes with K, but otherwise arbitrary, also 
is a Jf - odd. Let us remark for the further application, 
that the following useful relation holds in the framework 
of conditions of the above theorem 



E • 



1 



V xl-lxV 



(10) 



Now one can proceed to the second stage of our prob- 
lem - we wish to construct the K-odd operator A, which 
commutes with H . It is clear that there remains large 
freedom according to the above mentioned remark about 
O - operator - one can take O into account or ignore it. 

Wc have the following physically interesting vectors 
at hand which obey the requirements of our theorem. 
They are: 

r - unit radius vector and p - linear momentum vector. 

Both of them are perpendicular to I. Constraints of 
this theorem are also satisfied by Laplace-Runge-Lenz 
vector A = f - [p x I — I x p\, but this vector is 
associated to the Coulomb potential. Hence, we abstain 
from its consideration for now. We can remark that (E • 
A) is not an independent structure. It is expressible by 
two other structures, e.g. (E • i*) = (E • + :^l3K{'S-p) 
. Therefore, we choose the following K-odd terms: 



(E-r) and K{i:-p) 



(11) 



As it turns out inclusion of K into the second term 
of (11) is necessary for obtaining our final result. Let's 



remark that both operators in (11) arc diagonal matri- 
ces, while the Hamiltonian (1) is nondiagonal. Therefore, 
in commuting of (11) with nondiagonal terms appear as 
well. For instance, 



r 



(12) 



j:-r,H 

Therefore, we probe the following operator 

A = a;i (s • f) + ix2K (s • p) + ix^K-y^f (r) (13) 

Here the coefficients Xi {i = 1,2, 3) are chosen in such 
a way that A is Hermitian operator for arbitrary real 
numbers and / (r) is an arbitrary scalar function to be 
determined later. 

The commutator of A with H is calculated straight- 
forwardly, the result is 

[A, H] = xiffiKj^ + X2K(E-fj V (r) - 
- x^K (S • f) /' (r) - ix32mpK^^f (r) 

Equating the above expression to zero, i.e. requiring 
commutativity of our operator with the Dirac Hamilto- 
nian, we find 



K 



(E-f) [X2V' {r)-X3f{r)] + 
r T.I 1 

= 



(14) 



+2il3Kj^ y-mx3f{r) 



Here terms are grouped in a way that we have a diago- 
nal matrix in the first row, while the anti-diagonal matrix 
is in the second row. Therefore, the two equations follow: 



X2V' (r) = xsf (r) 



xamf (r) 



Xi 



(15) 



(16) 



Integration of the Eq. (15) with the requirement that 
functions / (r) and V (r) tend to zero when r — > oo, yields 



X3V (r) = Xsf (r) 
while the equation (16) gives 

/ (r) = 

X2 mr 



(17) 



(18) 



Substituting Eq. (18) into Eq (17) results in the fol- 
lowing potential 



V{r) = 



Xi 1 

X2 mr 



(19) 



Hence, in the very general framework we have shown 
that the only central potential for which the Dirac Hamil- 
tonian would have an additional symmetry (in the above 
mentioned sense) is a Coulomb potentiaL 
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Meanwhile, the relative signs of coefficients xi and X2 
may be arbitrary. Therefore we have a symmetry both 
for attraction and repulsion. If we take the Coulomb 
potential in the usual form 



r 

where a ~ Ze^ = Za , it follows 

1 

X2 = Xi 

ma 



(20) 



(21) 



In this case our symmetry operator (13) becomes 

A = xi{(j:-'r) - —K (^:■p) + —K-/^} (22) 
V / ma \ ) mr 

Number xi, as an unessential common factor may be 
omitted. Moreover, if we make transition to the usual 
Dirac a matrices according to the relation E = 7^a, then 
operator A can be reduced to the form 



A = 7^{a • r 



-Kj^ {H - Pm)} 



(23) 



which coincides precisely with the Johnson-Lippmann 
operator''. 

We mention here that if the potential in the Dirac 
Hamiltonian was a Lorentz-scalar (which means the 
change V — > PV ) then, while K still commutes with 
-ff , operator A does not commute anymore with H even 
for Coulomb potential. Thus, we are convinced that in 
this problem of supersymmetry, the Coulomb potential 
(as a fourth component of 4-vector, i.e. minimal gauge 
invariant switching) takes exceptional role - the super- 
charges and Witten algebra can be constructed only for 
this potential. What the real physical picture is stand- 
ing behind this? Remark that taking into account the 
Eq. (10) for V — p, one can recast our operator for the 
Coulomb potential in the following form 



A 



2ma 



p X I ~ I X p 



+ —Kf (24) 



One can see that in the non-relativistic limit, i.e. [3 
and 7^ — » , our operator reduces to 



Anr — 



2ma 



p X I — I X p 



(25) 



Note the Laplace- Runge-Lenz vector in the parenthesis of 
Eq. (25). Therefore, relativistic supercharge reduces to 
the projection of the Laplace-Runge-Lenz vector on the 
electron spin direction. Precisely this operator, Eq.(25) 
was used in the case of Pauli electroni. 

Therefore, we see that there is a deep relation between 
supersymmetry of the Dirac Hamiltonian and the sym- 
metry related to the Laplace-Runge-Lenz vector, which 
appeared already in classical mechanics and provides the 
closeness of celestial orbits. 

We can conclude that the hidden symmetry associ- 
ated to the Laplace-Runge-Lenz vector governs very wide 
range of physical phenomena from planetery motion to 
fine and hyperfine structure of atomic spectra. As for the 
Lamb shift, which is pure quantum field theory effect, its 
Hamiltonian, derived by radiative corrections to a photon 
propagator and photon-electron vertex function, does not 
commute with A operator and therefore spoils the above 
mentioned symmetry (supersymmetry). In other words, 
symmetry of the A operator controls the absence of the 
Lamb shift in the Dirac theory. 

In conclusion, one must also remember that the form 
of obtained symmetry operator is not unique. One 
can always replace A — *■ OA , where [0,H] = and 
[O, K] = . One can take, for example, O — f {K) - ar- 
bitrary regular matrix function of K . Moreover, SUSY 
in specific and mostly exotic models of Dirac equation 
(such as 2-1-1 dimensionsi^, non-minimal or anomalous 
magnetic moment couplingi^, squared equationi^*i^, 
etc.) are not excluded by our above consideration. 
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